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Abstract. We study the influence of many-particle interactions on a metal-insulator transition. We consider 
the two-interacting-particle problem for onsite interacting particles on a one-dimensional quasiperiodic 
chain, the so-called Aubry-Andre model. We show numerically by the decimation method and finite-size 
scaling that the interaction does not modify the critical parameters such as the transition point and the 
localization-length exponent. We compare our results to the case of finite density systems studied by means 
of the density-matrix renormalization scheme. 
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1 Introduction 



The metal-insulator transition (MIT) in disordered elec- 
tronic systems has been the subject of intense research 
activities over the last two decades and still continues to 
attract much attention. For free electrons in disordered 
systems [1] the scaling hypothesis of localization [2] can 
successfully predict many of the universal features of the 
MIT. However, the influence of many-particle interactions 
on the MIT is not equally well understood [3] and re- 
cent investigations of an apparent MIT in two-dimensional 
(2D) systems even question the main assumptions of the 
scaling hypothesis [4,5,6,7,8,9]. A simple theoretical ap- 
proach to the interplay of interactions and disorder is 
based on the two-interacting-particles (TIP) problem in 
ID random [10,11,12] or quasiperiodic potentials [13,14]. 
Furthermore, numerical results for spinless fermions at 
finite particle density have given additional insight [15, 
16,17]. In general, these investigations have shown that 
changes in the wave function interferences due to many- 
particle interactions [18,19] can lead to a rather large en- 
hancement of the localization lengths in ID and 2D [16, 
20,21]. 

The standard approach for computing localization len- 
gths in disordered, non-interacting systems is the transfer- 
matrix method [22] . It has been used for investigations of 
the enhancement of the TIP localization length in a ID 
random potential [12,23] where there is no MIT as all 
wave functions are always localized. Other numerical ap- 
proaches to the TIP problem have been based on the time 
evolution of wave packets [10,24], exact diagonalization 
[25] or Green function approaches [20,26,27]. 



In the single-particle case, the ID quasiperiodic Aubry- 
Andre model is known rigorously to exhibit an MIT for 
all states in the spectrum as a function of the quasiperi- 
odic potential strength /i [28]. The ground state wave 
function is extended for ^ < 1 and localized for /i > 1. 
The system at = 1 is critical: there the wave func- 
tions decrease algebraically, not exponentially as in the 
localized case. Recently, we examined this model for TIP 
by means of the transfer-matrix method together with a 
careful finite-size-scaling analysis [14] following earlier an- 
alytical work of Refs. [29,30]. We showed that the model 
for TIP exhibits an MIT as a function of /i at = 1 as 
in the single-particle case. Our finite-size-scaling results 
for onsite (Hubbard) interaction suggest that the criti- 
cal behavior, i.e., the value for the critical exponent u of 
the correlation length, is also not affected by the interac- 
tion [14]. However, it has been demonstrated [12,20] that 
a transfer-matrix-method approach applied to the TIP 
problem without finite-size scaling leads to unreliable lo- 
calization lengths, i.e., it systematically overestimates the 
TIP localization length A2 in finite-sized samples in the 
case of vanishing interaction (U — Q). In addition, simple 
extrapolations to infinite sample size [12,23] may lead to 
an underestimation of A2 [31]. An alternative approach, 
which does not suffer from the above problem, is based 
on the decimation method and has also been applied re- 
cently to TIP in a ID random potential [20]. This en- 
couraged us to reexamine the localization lengths for TIP 
in ID quasiperiodic potentials with Hubbard interaction 
with the decimation method. As we shall show in the fol- 
lowing, we find that the general conclusions of Ref. [14] 
remain valid, i.e., the MIT is not affected by the interac- 
tion. The critical properties of the single-particle transi- 
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tion at /ic = 1 are not altered within the accuracy of our 
calculation. One-parameter scaling is obeyed for onsite in- 
teraction strengths up to U = 10. 

As an independent extension of these low-density re- 
sults, Chaves and Satija [32] have studied a model of 
nearest-neighbor interacting spinless fermions [33] at fi- 
nite particle density in the same quasiperiodic potential 
by means of Lanczos diagonalization for small systems 
up to chain size M = 13. They have found evidence for 
a critical region. In order to reach much larger system 
sizes for interacting systems, one can employ the numer- 
ical density- matrix rcnormalization group (DMRG) [34]. 
With the DMRG the ground state properties in ID can 
be obtained very accurately [35,36]. In a recent paper [37], 
we studied the quasiperiodic model of Ref. [32] at various 
densities and interaction strengths V by DMRG. We com- 
pare the results with the present TIP data at the end of 
the paper. 

The paper is organized as follows. In section 2 we de- 
scribe the Hamiltonian of our TIP system and explain 
how to obtain the TIP localization lengths via decima- 
tion method. In section 3, we comment on the particular 
finite-size-scaling method employed and present the esti- 
mated critical parameters. We summarize and conclude in 
section 4. 



2 The TIP system and the numerical 
approach 

The Hamiltonian for TIP in the ID quasiperiodic potential 
of the Aubry- Andre model is given as 

H = \n,m){n +l,m\ + \n,m){n,m + 1\ + h. c. 

-I- |n, m) [/tin + ^J'm + U (n, m)] (n, m\ . (1) 

Here = 2ficos{a'm+(3) is the quasiperiodic potential of 
strength jj, with a/27r being an irrational number. /3 is an 
arbitrary phase shift and we choose a/2'K = {V^ — l)/2, 
i.e., the inverse of the golden mean. This value of a/27r 
may be approximated by the ratio of successive Fibonacci 
numbers — F„ = F„_2 + F^-i = 0, 1, 2, 3, 5, 8, 13, 
... — as is customary in the context of quasiperiodic sys- 
tems [38]. The Hubbard onsite interaction matrix U (n, m) 
is diagonal, i.e., U{n,m) = U5nm- The indices n and m 
correspond to the positions of each particle on a chain of 
length M. Now we use the decimation method [20,39] to 
construct an effective Hamiltonian for the diagonal of the 
M X M lattice along which the cigar-shaped TIP wave 
function has its largest extent [25,31]. The quantity of in- 
terest is the TIP localization length A2 defined by the TIP 
Green function G2{E) [26]: 

^ = -p^^ln|(l,l|G.lM,M)l. (2) 

For TIP in ID and 2D random potentials, this approach 
has led to high precision results supporting the proposed 



increase of the TIP localization lengths due to the re- 
pulsive interaction [20,21]. We remark that similar data 
have also been obtained for nearest-neighbor [26] and long- 
rangcd interactions [24]. 

The correlation length ^rx> for the infinite system may 
be obtained from the localization lengths A(M) for fi- 
nite system sizes by the using one-parameter scaling hy- 
pothesis Am — /(Af/?oo) [40] for the reduced locahzation 
lengths Am = \{M)/M. The MIT is characterized by a di- 
vergent correlation length ^00 (m) ^ IM~A*c1~'' [!]• In order 
to reliably extract the critical parameters from the calcu- 
lated values of X2{M) one may apply a finite-size-scaling 
procedure [22] that numerically minimizes deviations of 
the data from the common scaling curve /. The critical 
exponent v can then be extracted by fitting the ^00 ob- 
tained from finite-size scaling [41,42]. This method was 
used previously [14] for finding the critical parameters of 
the present model. 

Higher accuracy can be achieved by a method applied 
recently [41,42,43,44] to the MIT in the Anderson model 
of localization. We construct a family of fit functions which 
include corrections to scaling such as (i) nonlinearities of 
the dependence of the scaling variable on the quasiperiodic 
potential strength and (ii) an irrelevant scaling variable 
which accounts for a shift of the crossing point of the 
^m(a*) curves as a function of /i, i.e., 

^M = /(XrMi/-,XiM^) . (3) 

where Xi a-nd Xi ^I'c the relevant and irrelevant scaling 
variables, respectively. / is then Taylor expanded up to 
order m in terms of the second argument 

AM = Y,XiM^y~Ux.M''l , (4) 

71=0 

and each /„ is Taylor expanded up to order rir: 

fn = Y.^ni^rM"'' . (5) 

i=0 

Nonlinearities are taken into account by expanding Xt and 
Xi in terms oi u = {/ic — m)/Mc up to order rrir and mi, 
respectively, 

mr mi 
n=l n=0 

with 61 = Co = 1. The fit function is being adjusted to the 
data by choosing the orders ni, rir, mr, mi up to which the 
expansions are carried out. Of course, the orders have to 
be taken not too large to keep the number of fit parameters 
dni, bn, and Cn reasonably small. 

3 Numerical results for TIP 

We calculate A2 at energy E = for 20 values of the Hub- 
bard interaction, i.e., [/ = (the non-interacting single- 
particle case), 0.1, . . ., 0.9, 1, 2, . . ., 10 for 6 system sizes 



Andrzej Eilmes et al.: Localization properties of TIP in a quasi-periodic potential with a MIT 



3 



M = 13, 21, 34, 55, 89, 144. For [/ = and 1, we also 
have data for M = 233 and 377. The quasiperiodic po- 
tential strengths // is chosen close to the single-particle 
transition at /ic ~ 1 and ranges typically from 0.95 to 
1.05. As in Ref. [14] we average the results over differ- 
ent randomly chosen phase shifts /3 in order to reduce the 
fluctuations. The number of /3 values used in this averag- 
ing ranges from 5000 for M = 13 to 1000 for M = 377. In 
order to perform the non-linear fit necessary for the finite- 
size-scaling procedure as outlined in section 2, we use the 
Levenberg-Marquardt method [44, 45] . As the decimation- 
method data like the transfer-matrix-method results 
[14] — are still rather noisy we have to suitably limit the 
ranges of the quasiperiodic potential strength fx and/or 
the system sizes M used for fitting the data. 

For U = and 1 , which were examined by the transfer- 
matrix method in detail [14], the best fit is obtained for 
Tir = 3, rii = 2, nir = 2 and nii = 1. For U = we 
use the data for ^ ranging from 0.96 to 1.01 and M = 
55, 89, 144, 233, and 377; for [/ = 1 we use all system sizes 
M = 13, 377 and 0.97 < n < 1.05. Figures 1 and 
2 show the resulting TIP localization lengths for J7 = 
and 1. Also shown are the fits of the finite-size-scaling 
curves to the data as given by Eq. (3) for U = and 1, 
respectively. We find that for both U values, there is an 
apparent transition close to /ic = 1. For the case U = 0, 
we also observe a systematic shift of the crossing point 
with increasing system sizes necessitating the inclusion of 
an irrelevant scaling variable as discussed in section 2. 
The transition point is not so clearly distinguished for 
U = 1, albeit the different behavior for /i < 1 and /i > 1, 
namely the increase and decrease, respectively, of Am with 
increasing M, is clearly seen. 

The corresponding plots of the scaling curves are dis- 
played in the Figs. 3 and 4. The scaling curves are much 
better than reported previously [14] for the transfer-matrix- 
method data. The critical parameters can consequently be 
estimated to be /ic = 0.989 ± 0.001, u = 1.00 ± 0.15 for 
[/ = and /ic = 0.997 ±0.001, ly = 1.19 ±0.16 for U ^ 1. 
The irrelevant scaling exponents are close to y = 1.8 ± 0.2 
and y = 0.15 ± 0.1 for [7 = and 1, respectively. Note 
that the quoted errors correspond to the standard devia- 
tions estimated from the non-linear fit procedure. In this 
way the accuracy is significantly overestimated. Since it is 
apriori not clear, which values T?,i, n^-, m,-, nii to use. we esti- 
mate the true errors from a comparison of various fits with 
different ni, n^, m,^. mi. Even in the case of extremely high 
precision data close to the MIT in the Anderson model of 
localization, this has been shown [44] to increase the error 
by one order of magnitude. Therefore we conclude that 
the interaction strength U for TIP does not influence the 
MIT in the quasiperiodic potential within the accuracy of 
the present calculation. 

Further results for larger U values are collected in Ta- 
ble 1. The expansion orders n\,nr,mr,mi, the system sizes 
and ranges of the quasiperiodic potential strength have 
been chosen in order to minimize the statistics and 
to get the most convincing scaling fit. Furtlu^rmore, one 
has to check that various initial parameters (a„,6„,c„) 



converge to the same values of the critical quasiperiodic 
potential strength and the critical exponent v. Figures 
5 and 6 show the values obtained in this way. For almost 
all cases the critical quasiperiodic potential strength /ic 
remains close to 1, the only exceptions are U = Q and 
0.1, when ji^ = 0.99 and 0.98, respectively. However, since 
we know that the transition in the single-particle case is 
exactly at /Xc = 1 [28], this observation can be used to 
estimate the true error of the estimate for /ic. Thus com- 
paring with the /Xc estimates for [/ ^ 0, we find that the 
errors calculated within the non-linear fitting procedure 
are significantly underestimated as discussed above. We 
therefore conclude that within the accuracy of our cal- 
culation there is no change of the critical quasiperiodic 
potential strength /ic for the Hubbard interaction in the 
range < ?7 < 10. The same argument leads to the con- 
chision that within the error bars the critical exponent v 
does not change with the Hubbard interaction strength 
and is close to 1. This is an agreement with the previous 
results obtained by the transfer- matrix method and finite- 
size scaling [14] . We stress that the critical exponents can 
only be obtained with much less accuracy than the tran- 
sition point /ic as shown in Table 1. 



4 Conclusions 

In this work, we have studied the interplay of disorder 
and interactions for a quantum system at very low den- 
sity (TIP). We calculated the pair localization lengths 
for a quasiperiodic potential and Hubbard interaction by 
means of the decimation method and extracted the criti- 
cal parameters from the fit using the one-parameter scal- 
ing hypothesis. For both non-interacting particles as well 
as onsite interaction we obtain the value of the critical 
quasiperiodic potential strength /ic = 1 and the critical 
exponent s=; 1 in agreement with the previous results of 
transfer-matrix-method calculations and finite-size scaling 
[14]. The results for f/ > 1 show that this conclusion re- 
mains valid also for much stronger interactions. 

Let us briefly compare these results to the finite den- 
sity situation. For N interacting spinless fermions on a 
ID ring of circumference M with Aubry- Andre onsite po- 
tential /i and nearest-neighbour interaction V it is possi- 
ble to treat system lengths up to about M w 100 — 200 
using the DMRG. We applied [37] the finite lattice algo- 
rithm for non-reflectionsymmetric models as described in 
[46]. For a system of free fermions at finite density like 
p = 1/2 (incommensurate compared to the wave vector of 
the quasiperiodic potential — an irrational multiple of tt), 
we reproduced [37] the expected transition at /ic = 1 in 
agreement with Refs. [14,32]. For attractive and repulsive 
interactions at p = 1/2 the numerial results are available 
for only two system sizes (M = 34 and 144), therefore con- 
clusions about these regimes appear rather speculative. At 
commensurate densities pi « lim„^oo ^n-j/-Fri ~ 0.618, 
0.382, 0.236, and 0.146 corresponding to i = 1, . . ., 4 

and in the repulsive regime (nearest-neighbour interac- 
tion V > 0), the ground state is localized for /i > [37] 
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in agreement with previous studies for disordered and pe- 
riodically disturbed systems [15,36]. The above increase 
of the localization lengths as predicted by the arguments 
for TIP [10] is most likely too small [16] to bo detected 
by the present accuracy. For attractive interactions V, all 
densities pi and — > 0, the system shows a Peierls-Hke 
transition from insulating to metallic phase at V « — 1.4 
[37] in agreement with the weak-coupling renormalization 
group treatment [47] of spinless fermions on a Fibonacci 
lattice. 

In conclusion, we have studied the influence of inter- 
actions on an MIT in a quasipcriodic model in ID. Our 
results suggest that the delocalization found for low den- 
sity TIP in the localized phase cannot simply be extrap- 
olated to the finite-density situation. At finite densities, 
other effects such as a Peierls-like commensurability be- 
come important and dominate the transport properties. 

We thank M. Leadbeater for help with the decimation method 
and C. Schuster for stimulating discussions. We gratefully ac- 
knowledge the support of the SMWK and the Deutsche For- 
schungsgemeinschaft within Sonderforschungsbereich 393. 
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Table 1. Values of the critical quasiperiodic disorder strength /Xc and the critical exponent v obtained by the non- linear fit for 
various U values. The first row for each U gives values and the orders ni, mi, used in the expansion (4-6), for which the best 
fits have been obtained. In all cases we find nr = 3 and mr = 2. For // and M the range of the values which were used in the fit 
is given. The second row contains values of the critical parameters obtained from a weighted average of fits for various choices 

of Hi and mi. 
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Fig. 1. Reduced localization lengths Am versus quasiperiodic 
disorder strength fx for U = 0. For clarity, only error bars for 
M = 55 and 377 are given. The lines are the fits to the data 
given by Eq. (3). 
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Fig. 2. Reduced localization lengths Am versus quasiperiodic 
disorder strength n for U = 1. For clarity, only error bars for 
M — 377 are given. The lines are the fits to the data given by 
Eq. (3). 
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Fig. 3. Scaling function (solid line) and scaled data points for 
U = 0. For clarity only every 3rd data point is represented by 
a symbol. 
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Fig. 4. Scaling function (solid line) and scaled data points for 
U — 1. For clarity only every 3rd data point is represented by 
a symbol. 
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Fig. 5. The critical quasiperiodic potential strength /ic versus 
Hubbard interaction strength U. Error bars mark the errors 
resulting from the Levenberg-Marquardt method of the non- 
linear fit. 




Fig. 6. The critical exponent f versus Hubbard interaction 
strength U. Error bars mark the errors resulting from the 
Levenberg-Marquardt method of the non- linear fit. 



